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In modern treatises on Mathematical Physics it is customary to de- 
fine the derivative of a scalar function, taken at a given point in space 
in a given direction, in a manner which emphasizes the fact that this 
derivative is an invariant of a transformation of codrdinates. Accord- 
ing to this definition,! if through the point P a straight line be drawn 
in a fixed direction (s), if on this line a point P’ be taken near P so that 
PP’ bas the direction s, and if w,, up» be used to represent the values 
at these points of the scalar point function uw, then if the ratio 


Up—Up . 

PP (1) 

approaches a limit as P’ approaches P, this limit is called the derivative 

of uw, at P, in the direction s. If w happens to be defined in terms of a 

system of orthogonal Cartesian coérdinates, x, y, z, and has continuous 

derivatives with respect to these codrdinates everywhere within a 

certain region, the limit just mentioned exists in this region and its 
value is 


Ou 0 
008 (a; 8) + 5-008 8) 008 (2) 


1 Hamilton, Elements of the Theory of Quaternions; Tait, Elementary 
Treatise on Quaternions; Gibbs, Vector Analysis; Maxwell, Treatise on 
Electricity and Magnetism; Webster, Dynamics of Particles and of Rigid, 
Elastic, and Fluid Bodies; Jeans, Mathematical Theory of Electricity and 
Magnetism; Lamé, Lecgons sur les Coordonnées Curvilignes; Peirce, Theory 
of the Newtonian Potential Function; Generalized Space Differentiation of 
the Second Order; Czuber, Wienerberichte, 101,, 1417 (1892); Boussinesa, 
Cours d’ Analyse Infinitésimale; H. Weber, Die Partiellen Differential-Gleich- 
ungen der Mathematischen Physik. 


s 
ke, 
af 
# 
tes 
oe 
2. 
rf 
; 


340 PROCEEDINGS OF THE AMERICAN ACADEMY. 


Of all the numerical values which the derivative of w can have at a 
given point, the greatest is to be found by making s normal to the level 
surface of wu which passes through the point. This maximum value, 


du\? du\? Ou \? 

is usually regarded as the value at the point of a vector point function 
called the gradient vector of uw, the lines of which cut orthogonally the 
level surfaces of w, and the components of which parallel to the coérd- 
inate axes are 

Ou Ou Ou 

(4) 
This vector is, of course, lamellar. 

The value of the tensor of the gradient vector is often called simply 
the “gradient” of w and is denoted by h,. If at any point a straight 
Jine be drawn in the direction (”) normal to the level surface of wu, 
in the sense in which ~ increases, and if a length h, be laid off on this 


line, the projection, 
hy * COS (n, 8), (5) 


of this length on any other direction (s) is numerically equal to the 
derivative of « in the direction s. : 

Most physical quantities — such as temperature, barometric pres- 
sure, density, inductivity — present themselves to the investigator as 
single valued point functions, which, except perhaps at one or more 
given surfaces of discontinuity, are differentiable in the sense just 
considered. 

It is often desirable to differentiate a scalar function, w, at a point, 
in the direction in which another scalar function, v, increases fastest, 
and if (wu, v) represents the angle between the gradient vectors of w and 
v at the point, the derivative is evidently equal to 


hy * COS (u, v). (6) 


It frequently happens that in a question of maxima and minima, 
one wishes to determine the greatest (or the smallest) value which a 
quantity U7 may have, subject to the condition that another quantity 
V shall have a given value (V,). If these quantities can be represented 
by point functions, the problem geometrically considered requires one 
to find the parameter of a surface of the constant U7 family, which is 
tangent to the surface of the V family upon which V is everywhere 
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equal to V,; but at the point of tangency, the derivative of the function 
U in any direction in the tangent plane of the V surface is zero, that 
is, the normals to the U and V surfaces coincide, so that 


du Ou Ou 


dy 
(7) 


and these familiar equations usually furnish some general information 
about the problem independent of the value of V, As an extremely 
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Figure 1. 


simple example we may take the familiar problem concerning the rela- 
tive dimensions of an open tank of square base (a X 2) and height y, 
which shall hold a given quantity (V=.2*.y) of water and have the small- 
est wet surface (7 =2?+42y). Here we have the curve D of the V 
family, which has the given parameter, V,, and are required to find that 
member of the P, Q, #, S family which touches D. The equation (7) 
becomes in this case 2y=2, and it appears (Figure 1) that the curves 
of the two families which pass through any point of the line OW are at 
that point tangent to each other. 

It is sometimes necessary to differentiate a point function, wu, at a 
point P, in the direction of the line through the point, along which 


/ 
> 
MEE 
\ 
tary 
& 
she 
° 
s 


342 PROCEEDINGS OF THE AMERICAN ACADEMY. 


two other point functions, v, w, are constant ; that is, along the line 
=v,w=ew, If 


av ax av av av av 
Oy dz dz Ox Ox Oy 
L= », M= , N= (8) 
aw dw aw 
dy dz dz Ox Ox Oy 


and if Rk? = L? + M* + N*—which is equal to h,? - hy’, if v and w 
are orthogonal — this direction is defined by the cosines L/R, M//R, 
N/R, and the derivative required is 


1 Ou Ou Ou 
(9) 


If the maxima and minima of the function uw = / (2, y, z) are to be 
found under the condition that the functions v, w shall have given 
numerical values, the derivative of w taken in the direction in which v 
and w are constant must be made to vanish. Thus, if 


2, 
and if the conditions are 
ayz=c and 
equation (9) yields immediately the required relation 
(zy + = 0. 


When /’ (wz) is positive, the direction of the gradient vector of / (w) 
coincides with that of the gradient vector of w itself: these directions 
are opposed when ,/’ (w) is negative. The tensors of both vectors are 
always positive. If 


w=f(u), (u)}* - and cos (w, s) = cos (u, s): 
in particular, when 


= 1/u, hy =h,,/u? and cos (w, s) = — cos (u, 8), 
so that 
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If w is the distance (r) to a point on a curve (s) from a fixed point 
outside the curve, 


Any function of the complex variable (az + by + izWa? + 6?) has a 
gradient identically equal to zero, but every differentiable real point 
function has a gradient in general different from zero. ‘The gradient of 
a function may be constant throughout a region of space: if the 
gradient of w is constant, the surfaces upon each of which w is constant 
form a parallel system. If the gradient of a function, w, is either con- 
stant or expressible in terms of w, any differentiable function of w has 
a gradient either constant or expressible in terms of w. If the gradient 
of wu is expressible in terms of w alone [h,, = /(w) |, it is possible to form 


a function, a f FO , of uw the gradient of which shall be constant. If 


h,, is neither constant nor expressible in terms of w, no function of 
exists the gradient of which is expressible in terms of wv. The functions 
2), v= sin + 2y — 3z), w=sin (5a — 4y — z) 
illustrate the fact that the gradient of each of three orthogonal point 
functions may be expressible in terms of the function itself. 

If the gradient of each of two orthogonal point functions, u, v, were 
expressible as the product of a function of w and a function of ¥, so that 
h, = U,- Vi, and hy = U2- V2, it would be possible to form two func- 


tions LF | of w alone and of alone, respectively, the 


gradient + ‘als of which would be expressible in terms of the other. 
If the gradient vectors of two functions have the same direction at 
every point of space, one of these functions is expressible in terms of 
the other. If the gradients of two real functions, wu, v, are everywhere 
equal while the directions of their gradient vectors are different, 


r) 


and the functions [w + 7], [w— v] are orthogonal, as are + v), 
(u — v), where and / are any differentiable functions. If wand v 
are orthogonal functions, the functions + /(v)], —/(@)] 
have gradients numerically equal to each other at every point. 

T'wo scalar point functions, the level surfaces of which are neither 
coincident nor orthogonal, may have gradients each of which is ex- 
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pressible in terms of the other: the gradient of v = $2*— 4 ry? is 
equal at every point of the zy plane to the square of the gradient of 
u=2*—y’. If u and v are orthogonal functions of z and y, the 
product of their gradients is equal to the Jacobian, 


Ou Ou Ov. 
Ox Oy dy Ox 
The differential equation 


du \? du \? du \? 
(az) + Gy) + (Gz) 
which leads to systems of parallel surfaces, is of standard form. Its 
complete integral is 


where a, b, d are arbitrary constants, and from this the general 
integral may be obtained in the usual manner. 


If a direction s be determined at every point of a given region, 7, 
by some law, the derivative of the function « becomes itself a scalar 
point function in 7’, and if this is differentiable, it may be differentiated 
at any point in any direction, say s. It is usually convenient to 
define s by means of three scalar point functions, 4, m, m, the sum of 
the squares of which is identically equal to unity, and which represent 
the direction cosines of s. In this connection it is well to notice that 
if s has the direction at P of the tangent of a continuous curve which 
passes through the point, if P’ be a point near P on the tangent and 
P” a point near P on the curve, and if U is any differentiable scalar 
point function, 

pu — Up Up — Up 
re” rr’ 


have the same limit, as P’ and P” approach P, that which has been 
defined as the derivative of U at P in the direction s. If, then, 
du/ds is differentiable 


(Ou Ou Ou Ou 
Ou Ol du Om du On 
dz da’ 


(11) 
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and 

(12) 


If s’ is a direction defined by the cosines /’, m’, n’, 


~~ u u 
Ouf,, al , Om 
Ouf. On 


and it is clear that the order of differentiation is usually not com- 
mutative. Derivatives of this kind are often found in differential 
equations of orders higher than the first which define functions in 
terms of simple curvilinear coérdinates. 

If for instance spherical codérdinates are to be used, the second 
derivative of w taken in the direction in which @ increases fastest is 


ru 2 
29 cos*d + — 29 24 
COSY cos -cos’6 sin? y + sin -cos? sin cos 
2 07u 20°u 
— sind cos 6 cos  — cosGsin  — - sin cos ¢ 


Oy Oz 


and this, which contains derivatives of the first order, is in sharp con- 
trast to the second derivative of w taken in the direction 7, which is, 


20 
+ “sin @ cos sing + 


by de - sin 8 cos 6 cos ¢. (15) 
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Sometimes s and s’ are fixed directions so that J, m, n, I’, m’, n’, are 
constants throughout 7) and in this case the coefficients of du/dz, 
du/dy, du/dz in (12) and (13) vanish. The mutual potential energy W, 
of two magnetic elements, W/, M’, of moments, m, m’, can be written 


in the form 
— 
m-m >) (16) 


where 7 is the distance 1/M’ and s, s’ are the directions of the axes 
of the elements. The force (due to the second magnet) which tends 
to move the first magnet in the direction of its own axis is then 


and these differentiations assume that the direction cosines of s and 


s’ are constants. 
In general, if s is the direction perpendicular to the level surface of 


u, and if / is the scalar point function which gives the value of du/ds, 


07u Oh Ou . Oh Ou . Oh Ou 


In the case of oblique Cartesian codrdinates in a plane, 2 increases 
fastest in a direction which is not perpendicular to the line along 
which it is constant. Ifthe angle between the codrdinate axes is », 


Ou Ou Ou 

h,-cos (a, h,), h,-cos (y, hu), h,,- cos (s, h,), 
Ou _ Ou sin s) du sin (2, 8) (19) 
és sinw dy sinw 


It is frequently necessary to differentiate one point function, U, with 
respect to another, y, and the process usually appears in the form of a 
kind of partial differentiation. If, for instance, U is to satisfy a differ- 
ential equation in terms of a set of orthogonal curvilinear codrdinates 
of which w is one, the derivatives of 7 with respect to w are to be taken 
on the assumption that the other codrdinates remain constant. This 
large subject bas been treated exhaustively in the many works on 


a 
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orthogonal coérdinates which have been published since Lamé’s clas- 
sical treatise 2 appeared. 

Given a function, w, it is, however, not generally possible to find a 
system of orthogonal functions of which w shall be one, and it is often 
convenient for a physicist to differentiate a physical function, U, with 
respect to another, w, without considering the existence of any other 
related functions. A physical point function has a value at every point 
in space which is not altered by changing the system of codrdinates 
which fix the position of the point, and it is well to define the deriva- 
tive of UY with regard to wu in a manner which shall emphasize the fact 
that the derivative is an invariant of a change of codrdinates and which 
shall not assume that two functions (v, w) can be found orthogonal to 
each other and to w. When U and w are considered by themselves and 
not regarded as codrdinated of necessity with other similar quantities, 
it is usually, if not always, the case that a “normal” derivative? is 
required. 

The normal derivative, at any point, P, of the differentiable scalar 
point function U, with respect to the differentiable scalar point function 
u, may be defined as the limit, when PP’ approaches zero, of the ratio 


Up — Up 
upr— up’ 


(20) 


where /” is a point so chosen on the normal at P of the surface of con- 
stant w which passes through P, that wp-— up shall be positive. If 
(U, u) denotes the angle between the directions in which U and w in- 
crease most rapidly, the normal derivatives of U with respect to w and 
of w with respect to U may be written 


(U,u) and (U, u). (21) 
h,, hu 


If hyo= hy, these derivatives are equal. An example of this is the 
equality of dn/dr and dr/dn in a familiar application of Green’s 
Theorem, where » and 7 represent the normal distance from a given 
surface and the distance from a given fixed point respectively. If UV 
and w happen to be expressed in terms of a set (2, y, z) of orthogonal 


2 Lamé, Lecgons sur les Coordonnées Curvilignes et leur Diverses Appli- 
cations; Salvert, Mémoire sur l’Emploi des Coordonnées Curvilignes; Dar- 
boux, Legons sur les Systémes Orthogonaux et les Coordonnées Curvilignes; 
Goursat, Cours d’Analyse Mathématique. 

8 Peirce, Short Table of Integrals, Theory of the Newtonian Potential 
Function; Generalized Space Differentiation of the Second Order. 
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Cartesian codrdinates, the normal derivative of U with respect to wu 
can be written 
aU du, du, du 
Ox Ox Oz (22) 


and it is easy to see that this is equal to the ratio of the derivatives of 
U and u taken in the direction in which w increases most rapidly. 

It is occasionally instructive to use the conception of normal differ- 
entiation in studying some of the general equations of Physics: thus 
in the uncharged dielectric about an electric distribution, the potential 
function, V, is connected with the inductivity of the medium, pu, by the 
familiar equation 


in which » is to be regarded as a point function discontinuous in gen- 
eral at each of a given ‘set of surfaces at every point of which an equa- 
tion of the form 


+ = 0 (24) 
1 2 


is satisfied. Now (23) may be put into the form 


d lo 
av 


and according to Lamé’s condition, the second term is a function of V 
only, if the level surfaces of V are possible level surfaces of a harmonic 
function. 

It is easy to make from (25), by inspection, such simple deductions 
as those which follow in this paragraph. If V is harmonic, either the 
dielectric is made up of homogeneous portions separated from one an- 
other by equipotential surfaces, or the level surfaces of u and of V are 
everywhere perpendicular to each other. If V, though not harmonic, 
satisfies Lamé’s condition [V?(V)/ h,?= F'(V )] the level surfaces of the 
inductivity are equipotential ; and if the level surfaces of V and p» are 
identical, V satisfies Lamé’s condition. If when the plates of a con- 
denser are kept at given potentials, the level surfaces of the inductivity 
of the dielectric are equipotential, the value of the potential function in 
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the dielectric would be unchanged if » were changed to 2.y, where Q is 
any scalar point function orthogonal to V. If the continuous dielectric 
of a condenser in which the level surfaces of the inductivity, », are 
equipotential be changed so as to make the new potential function 
between the plates a function [V’ =/(V)] of the old, the new induc- 
tivity must satisfy an equation of the form p’ =OQ.u//’(V). If the V 
and the » surfaces are neither coincident nor orthogonal, V cannot be 
harmonic, and if V is given and one value of the inductivity found, no 
other value of the inductivity with the same level surfaces as this can 
be found except by altering the old value at every point in a constant 
ratio. If V does not satisfy Lamé’s condition, a new value of the 
inductivity found by multiplying the old value by any point function 
orthogonal to V, will yield the same value of V, but the level surfaces 
of the inductivity will be altered. Ifthe V and the » surfaces are not 
coincident, no change of the inductivity which leaves its surfaces un- 
changed can make these surfaces equipotential. 

If a mass of fluid, the characteristic equation of which is of the form 

=f (p, 7’), is at rest under the action of a conservative field of force 
the components of which are X, Y, Z, 


a 
~=p-Y, 


op 
0z 
It follows immediately from these equations that p and V must be 
colevel, and the normal derivative of p with respect to V shows that 
equilibrium is impossible unless the distribution of temperature is such 
that the equipotential surfaces are also isothermal. 

If the scalar point function, W, is expressed in terms of the three 
orthogonal point functions, uw, v, w, the square of the gradient of W is 
well known to be equal to 


, (aw? aw\? (aw? 


If the vector point function Q is expressed in terms of u, v, w, the 


divergence of @ is equal to 


If the normal derivatives of « and v with respect to w be denoted by 
D,,u and D,,v, it follows from the definition that 
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Dy (U + v”) = Dyu + Du" = n-u"—- 


— 


(u-v) = v- Dou + u- Dyn, 
Dw Jf (u) = Dw (u). 


The normal derivative of w with respect to v is a scalar function 
which, if differentiable, has a normal derivative with respect to v, and 


since by definition 


D «= (27) 


_ 1 ( Ov , Ohy OV , Dhy Dv 


we may write 


T ae- Oy Ox dy dy-dz > Oz" Oz Ox 
1 dv dh, dh, dv Oh, 


dh, dv Oh, 


1 Ow Ov . Ow Ov dw dav 
= | bat aa t ay dy dz 


+ 


hy” -hy* | Ox-dy\ dy Ox Oy], Oy dy dz 
(dw dv . dw dw 
1 duf av dw av Ow ov dw 


h, ax da\ da dy ay er 
dy? dy ay Oy-Ox dy 
hy dy dy\dy dy dz dz" Ox Aa 
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hy? + hyp? dz Ou =Az-dy dy 
2 ou Ow dhy dw dh, (30) 


O2 Oz\ Oz Oz Ox Ox Oy dy 


It is evident that D,D,,u is usually quite different from D,D,u. 


In the transformation of a partial differential equation from one set of 
independent variables to another set which does not form an orthogonal 
system, derivatives occur which are not normal in the sense of the last 
paragraphs. Ifa mass of fluid is in motion under the action of given 
forces, it is usually convenient either to express the orthogonal codrdi- 
nates of a particle which at the time ¢ has the position (2, y, z) in terms 
of ¢ and the codrdinates 2,, ¥, %, which the same particle had at the 
origin of time, or to express 2, Y, 2%, a8 functions of z, y, 2, t. 


X =Si (2, Y; t), Yo = fo (2, Y,%,t), %=S (a, y, % t). (31) 


In this case, it frequently happens that the level surfaces of A, A, /s, 
are not orthogonal. According as we use the “historical” or the “sta- 
tistical’ method of studying the motion, we shall express the pressure 
and the density in terms of 2, Y, %, ¢, or in terms of z, y, z, ¢. Sup- 
pose the second method to have been chosen, and dp/dzr to have 
been found by the aid of Euler’s Equations of Motion and the Equation 
of Continuity, and suppose that dp / 02, is needed. We shall then have 


dp _ Op dw | Op dy | Op a (32) 
oy OX, 0z OX, 


If with the help of (31) we find the values of the determinants 


Oy Oz (102 Ox Ox Oy 
L= , N= (33) 
Oy Oz dz Ox Ox Oy 


and put 


P=27+ M+  N* 
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we may write the results of differentiating all the equations of (31) 
with respect to x, Yo, %, in the form 


so that 


(35) 


and this is evidently equal to (9), the ratio of the directional deriva- 
tives of p and 2, taken in the direction (s) in the (2, y, 2) space in which 
both y, and z, are constant. If (s, py), (s, 2) represent the angles between 
s and the directions of the gradient vectors of p and a respectively, 


op _ hy * COS (s, Pp) (36) 


02, cos (8, 2)’ 


It is convenient, therefore, to define the derivative of a scalar point 
function, «, with respect to another scalar point function, v, at any 
given point in any direction (s), as the ratio of the directional deriva- 
tives of uw and v taken at the point in the direction s. 

Derivatives of this kind which frequently appear in two dimensional 
problems in Thermodynamics and in Hydrokinematics, usually involve, 
as has been said, a transformation from one set of codrdinates to an- 
other which is not orthogonal. 


JEFFERSON PuysicaL LABORATORY, 
CAMBRIDGE, Mass. 
December, 1909. 
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